Supersymmetry identifies molecular Stark states 
whose eigenproperties can be obtained analytically 
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We made use of supersymmetric (SUSY) quantum mechanics to find a condition under which 
the Stark effect problem for a polar and polarizable closed-shell diatomic molecule subject to 
collinear electrostatic and nonresonant radiative fields becomes exactly solvable. The condition, 

2 

Acj = 4( - J7 ^ +1 - ) 2 > connects values of the dimensionless parameters uj and Alv that characterize the 
strengths of the permanent and induced dipole interactions of the molecule with the respective 
fields. The exact solutions are obtained for the |J = m,m;uj, Aw) family of "stretched" states. The 
field-free and strong-field limits of the combined-fields problem were found to exhibit supersymme- 
try and shape-invariance, which is indeed the reason why they are analytically solvable. By making 
use of the analytic form of the |J = m,m;u),Auj) wavef unctions, we obtained simple formulae for 
the expectation values of the space-fixed electric dipole moment, the alignment cosine, the angular 
momentum squared, and derived a "sum rule" which combines the above expectation values into a 
formula for the eigenenergy. The analytic expressions for the characteristics of the strongly oriented 
and aligned states provide a direct access to the values of the interaction parameters required for 
creating such states in the laboratory. 

PACS numbers: 32.60.-H, 33.90.+h, 33.15.Kr,11.30.Pb, 37.10.Vz, 03.65.-w, 03.65. Ge 
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I. INTRODUCTION 

Whether a problem in quantum mechanics is exactly solvable is closely related to its supersymmetry (SUSY) [HE]- 
The foundations of supersymmetric quantum mechanics have been worked out by Witten in 1981 as an example of 
SUSY in zero-dimensional field theory [3J. Soon thereafter, supersymmetric quantum mechanics rapidly evolved into 
a new branch of mathematical physics [TJ HI 0] and reached a highpoint in 1983 - then unnoticed - when Genden- 
shtein established a connection between supersymmetry and exact solvability [5]: he demonstrated that Schrodinger's 
equation is exactly solvable if the potential and its superpartner exhibit shape-invariance. Whereas a supersymmetric 
Hamiltonian can be constructed for any potential whose ground-state wavefunction is analytic, shape invariance only 
exists for supersymmetric potentials that are interconvertible by a change of a parameter other than the integration 
variable itself [6 . Herein, we make use of the methods of supersymmetric quantum mechanics to arrive at exact 
wavefunctions and other eigenproperties of linear molecules subject to nonresonant electric fields in closed form. 

In our previous work on the molecular Stark effect, we showed that for polar molecules, combined collinear electric 
and nonresonant radiative fields can synergetically produce spatially oriented pendular states, in which the molecular 
axis librates over a limited angular range about the common field direction [7]. These directional states comprise 
coherent superpositions or hybrids of the field-free rotational states \J,m), with a range of J values but a fixed 
value of to, which remains a good quantum number by virtue of the azimuthal symmetry about the fields. This 
has proved an effective and versatile means to produce oriented molecules for applications ranging from molecule 
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optics and spectroscopy to chemistry and surface science [8HTU]. However, the eigenproperties of the Stark states in 
question had to be found numerically typically by diagonalizing a truncated Hamiltonian matrix. Here we show that 
supersymmetric factorization of the Hamiltonian yields exact wavefunctions \ J = m,m; w, Aw) in closed form for a 
particular ratio of the parameters w and Aw that determine the interaction strengths of the molecules with the static 
and radiative fields, respectively. 

The paper is organized as follows. In Section [nj we first introduce the molecular Stark effect problem for the case 
of collinear electrostatic and radiative fields. In Section |TO| we use methods of supersymmetric quantum mechanics to 
find solutions to the problem of molecules in combined fields in closed form. We derive a relation between strengths of 
the electrostatic and laser fields, w and Aw, at which the wavefunctions of the | J = m, m; w, Aw) states take a simple 
analytic form and show that SUSY furnishes isospectral partner potentials that can be realized with combined fields. 
In Sect ion [TV| we discuss shape-invariance of the supersymmetric partner potentials obtained and the general conditions 
for the exact solvability of the Stark effect problem. In Section [V] we investigate the field-free and strong-field limits, 
where the field strengths w and Aw approach, respectively, zero or infinity. We show that in both limits, the problem 
exhibits supersymmetry and shape-invariance and, therefore, can be solved exactly. Section |Vlj gives a summary of the 
exact closed-form expressions for the properties of molecules in fields obtained from the exact wavefunctions. These 
properties include the eigenenergy, the space-fixed dipole moment, the alignment cosine, and the expectation value of 
angular momentum. The availability of these otherwise hard-to-come-by properties in closed form allows to reverse- 
engineer the problem of finding the interaction parameters required for creating quantum states with preordained 
characteristics. The main conclusions of the work are summarized in Section [VII| An appendix surveys the concepts 
pertinent to exact solvability. 



II. MOLECULES IN COLLINEAR ELECTRIC AND RADIATIVE FIELDS 

We consider a 1 !) linear molecule with a rotational constant B, a permanent dipole moment \i along the internuclear 
axis, and polarizability components an and a± parallel and perpendicular to the internuclear axis. The molecule is 
subjected to an electrostatic field e combined with a nonresonant laser field of intensity /, whose linear polarization 
is collinear with e. With energy expressed in terms of B, the Hamiltonian takes the dimensionless form [llj : 



H = J 2 + V^ a (9), (1) 



with the angular momentum operator 



t2 1 d ( . „ d \ 1 d 2 

J = -sin^M Sm ^J-sT^W (2) 



and the interaction potential 



V^ a (6) = -uj cos 6- (Aw cos 2 6> + w ± ) (3) 

The dimensionless interaction parameters are given as w = \xejB and Aw = W|| — wj_, with uj\\.± = 27raii j_ I /{Be). 

The common direction of the collinear electrostatic and linearly polarized radiative fields defines an axis of cylindrical 
symmetry, chosen to be the space-fixed axis Z. The projection, to, of the angular momentum J on Z is then a good 
quantum number while J is not. However, one can use the value of J of the field-free rotational state, Yj tm (8, cf>), 
that adiabatically correlates with the hybrid state as a label, designated by J, so that —> Yj, m (6,(j)) for 

w, Aw — > 0. For arbitrary interaction strengths, the solution to the the Schrodinger equation with Hamiltonian ([T]) is 
an infinite coherent superpositions of the field-free wavefunctions, 

| J, to;w, Aw) = ]T c^(w»y Jmi (4) 
J 

whose expansion coefficients Cjj^ l (w, Aw) can be obtained by truncating series Q at some maximum value of J and 
then diagonalizing Hamiltonian |T]) in the finite basis set of the field-free wavefunctions. 

The axial symmetry of the problem allows to separate angular variables and express the dependence on the azimuthal 
angle 4> via the good quantum number to. In this case the Schrodinger equation with Hamiltonian can be recast 
as 



1 d ( . d 
suit 



sin 9d6 \ d9 



m = mo), (5) 
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where the three-dimensional effective potential is given by 
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K (3 a D) (°) = -ujcosO-Alo cos 2 9 (6) 

sin 

We note that all rotational levels are uniformly shifted by oj±. In what follows we use E = E^^ + uj± instead of the 
'true' molecular energy, £" MiQ . Moreover, since the Stark effect does not depend on the sign of to, we can define the 
projection of the angular momentum on Z as a positive quantity, to = \m\. 
By means of the substitution [T2"] . 

^(0) = /(0)(sin0)-i, (7) 
Schrodinger equation ([5| can be transformed to a one-dimensional form, 

d 2 



d6 2 f W 



/(*) = 0, (8) 



where 

^ ( ,a D) (0) = - w cos - A W cos 2 6—\—E (9) 

sin ' 



is a one-dimensional effective potential, which will be shown to play the role of one of the requisite superpartner 
potentials leading to ground state energy E — Eq. 

III. SUPERSYMMETRIC FACTORIZATION OF THE COMBINED-FIELDS HAMILTONIAN 

Here we invoke supersymmetry to find analytic solutions to eq. ([8]) and subsequently to eq. ([5]). Supersymmetry 
makes use of the first-order differential operators, A ± = +W(8), with W(9) the superpotential. The superpartner 
Hamiltonians are defined by 

H T =A ± A^ = -^+V^ m \6l (10) 

with the one-dimensional partner potentials V± D \0) = W 2 {9) ± W'(9). The superpartner Hamiltonians have the 
same energy spectra except for the ground state, and, if the eigenfunctions of one of the partner Hamiltonians Hzp 
are known, the eigenfunctions of the other can be obtained analytically via the intertwining relations [2 1^] . 
For the case of a molecule in combined fields we start from the following Ansatz for the superpotential: 

W{9) = acot(9) +q(9), (11) 

where the constant a and the function q(9) are to be determined. The first term in eq. (11) is the superpotential for 
a field-free rigid rotor, which is a special case of the Rosen-Morse I potential [2J. The superpotential W(9) yields the 

supersymmetric partner potentials V± D ^ (9): 

vi 1D) {9) = W 2 {8)±W'{9) = Q(Q y ^ + q 2 {9) ± q'{9) + 2aq{9) cot 9 - a 2 . (12) 
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By identifying the effective potential (|9| with V_'{Q) and substituting it into eq. (12), we obtain 



4 - u cos 6- Alu cos 2 9- - - E = +q 2 {9) - q'(9) + 2aq{8) cot 9 - a 2 , (13) 



1 _ _ o(o + l) 



sin 2 9 4 sin 2 (9 

which is satisfied for 



q(0) =/3sm9, (14) 
a = -(to +1/2), (15) 
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E = m(m + 1) - /3 2 , 



(16) 



and 



Alo = 



4(m + I)'- 



Substitution from eqs. (14) and (15) into eq. (11) leads to a superpotential 



W{9) = - ( m+- ) cot6> + /3sin6>, 



corresponding to the following pair of superpartner potentials: 



T/i 1D) (9) = m 4 - 2(3(m + 1) cos 6 - /3 2 cos 2 9 + f3 2 - m(m + 1) 
sin 9 



(17) 



(18) 



(19) 



Vi m) (9) 



(m + l) 2 - j „ 9 9 , 1 

4 - 2/3m cos 0-/3 2 cos 2 9 + f3 2 -m(m+ 1) - - 



sin 2 (9 



(20) 



We note that by identifying the effective potential ^ with V^ 1D ' I (#) instead, the result would have been the same for 
a = (m- 1/2). 

Thus, for each value of the good quantum number m there exists a pair of the supersymmetric partner Hamiltonians, 
ii_ and H + . which are isospectral and fulfill the following relations: 



which imply: 



H + (Af-(d))=E-(Af-(9)); 
H^f-(9))=E+(A^f-(9)) 



E n ~ E n+\i E = 



(21) 
(22) 

(23) 



f+(9) = (E- +1 )-^ 2 Af- +1 (9): f- +1 (9) = (£+)" 1/2 A* ft (9), 



(24) 



Afo(8) = 



The ground state wavefunction / (9) can be obtained from the superpotential ( 18 1 



/ -(#)=iVexp 



W(x)dx 



which yields the wavefunction in closed form, 

fo(9) = 7V(sin(9) (m+1/2) e Pcost 



Its normalization constant is given by 



N = 



T(m + 3/2) 



27T 3 /2r(m + l) Fi{] m + 3/2; /3 2 ) 



1/2 



(25) 



(26) 



(27) 



(28) 



where qFi(; a; z) is the confluent hypergeometric function [14 . Since the ground state wavefunction (27) is normaliz- 
able and obeys the annihilation condition (25 1, the supersymmetry obtained is unbroken [21 1151 US] . The top panels of 

Fig. 1 show the ID supersymmetric partner potentials V^ 1D ' 1 (6*) + Eq, the superpotential W{9), and the ground-state 
wavefunction fo(9) for \m\ = 0, 1, 2, 3 and different values of the interaction parameter (3. As f3 increases, the nodeless 
wavefunction fo(0) becomes strongly confined near the potential minimum. The potentials (19) and (20) come close 
to each other with growing m and eventually coincide in the limit of m 3> 1. 




FIG. 1: Top panels: ID supersymmetric partner potentials V_ (6) + Eq [black solid line], Vl (6) + Eo [red solid line], 
superpotential W(6) [green dashed line], and ground-state wavefunction fo(0) [blue solid line] for |m| = 0,1,2,3 and different 



values of /3. Bottom panels: 3D supersymmetric partner potentials, V^_ 3D \9) [black solid line] and Vj: J1J; (#) [red solid line] 
corresponding to the above ID potentials, along with the ground-state wavefunction of a molecule in combined fields, ipo(9) 
[blue solid line]. The eigenenergy of the ground state, Eq — m(m + 1) — /3 2 (in units of B), is shown by the black dotted line. 



r(3D) i 



By comparing the ID potentials V+ 1B ^ and V^_ 1B \ eqs. (19), (20 1, with eqs. (|5| and flij), we obtain the 3D potentials 
corresponding to a molecule interacting with the combined fields , 



2(m + l)/3cos0- /3 2 cos 2 



(29) 



V 



(3D) {0) = {m + lf _ ^ Q _ p2 cog2 ( 
sin 9 



(30) 



The ground-state wavefunction of the V^ 30 " 1 potential (corresponding to energy Eq — m(m + 1) — /3 2 ) then becomes 

4> {9) = N{-l) m \s\n6) m e ^ cose . (31) 

The phase factor (—1)™ leads to the correct asymptotic behavior of the wavefunction which, for (3 = 0, reduces to 
the ground state wavefunction of a rigid rotor with J — m, Y mtm (0,0). 

Eq. (30 1 coincides with the potential for a rigid rotor in the combined fields whose projection quantum number is 
m + 1 and whose interaction strengths are related by 



4m 2 



/3 2 



(32) 



Hence we have shown that, given a value of j3, the Hamiltonian of a molecule with a projection m of the angular 
momentum on the combined fields whose interaction parameters are related by eq. ( 17 1 has the same set of eigenvalues 
as a molecule with a projection m + 1 on the combined fields whose interaction parameters are related by eq. (32 1. 
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The bottom panels of Fig. [ljshow the 3D superpartner potentials, V^ 3D "*(0), for a molecule in the combined helds 
along with the ground-state wavef unctions, ipo(9), for \m\ — 0, 1. The permanent dipole term (oc cos#) is asymmetric 
with respect to 9 = ir/2 whereas the induced-dipole term (oc cos 2 9) is symmetric about 8 = ir/2 and gives rise to 
a double- well for an > a±. At large /3, w -C Aw and the coupling between the double- well's tunneling doublets 
results in localizing the ground-state wavefunction in the forward well of either the V^ 1D ^ or l/( 3£) ) potential, which 
for w -C Aw is mainly due to the polarizability interaction. At small /3, u> is comparable with Aw and the effective 
potential becomes skewed as the forward well becomes deeper at the expense of the backward well. For any j3 > 0, 
the molecule exhibits orientation (as opposed to alignment), which rapidly increases with f3. 

We note that for to = the ground-state level always lies deeper than the potential minima at 9 = tt, since 
V{9 = tt) — Eq = 2/3 > 0. Therefore, for any value of /3 the molecule remains confined to the potential minimum at 
9 = and exhibits orientation. For the "stretched" states of highly rotationally excited molecules, J = m S> 1, the 



superymmetric partner potentials ( 29 ) and ( 30 ) coincide 



IV. SHAPE INVARIANCE AND EXACT SOLVABILITY OF THE MOLECULAR STARK EFFECT 

PROBLEM 

Gendenshtein [5] demonstrated that the Schrodinger equation for any of the superpartner potentials, V-(9, a) and 
V+(#,a), admitting normalizable solutions, is exactly solvable if the potentials are translationally shape- invariant, i.e. 

V+(6,a ) + g(a ) = V-(6,a 1 )+g(a 1 ), (33) 

where the values ao and ai of the a parameter pertain to the n-th eigenstate with n = and n = 1, and the function 
g{a) is independent of 9. The parameter a\ is a function of ao, i.e. a\ = p(ao). 



If the superpartner potentials satisfy relation ( 33 1 , the eigenenergies and wavefunctions of ii_ can be obtained 
from m [17|: 



E 



n =g(p (n) M)~g(ao), (34) 



f-(6,a ) = At(0, a „_i)/ o -(0,a«), (35) 

where p^ n \ao) designates the function p(a) applied n times. 

Recently, Gangopadhyaya and Mallow demonstrated that condition (33) is equivalent to the following partial 
differential equation for the superpotential [TTJ HH] : 



2W(9,a) 



f dW(9,a) \ n dW(9,a) dg(a) 
' a) {—da— ) 2 ^9~ + ^a-=° (36) 



Unfortunately, equation ( 36 ) is not satisfied for the superpotential W(9, a) of eq. ( 18 ), neither for a = /?, nor for a = to, 

and, therefore, the partner potentials V_ D \9, a) and V^ 10 ^, a) are not shape invariant. Although, shape invariance 
is a sufficient, but not necessary condition for exact solvability |19j . Scrodinger equation for a molecule in combined 
fields is known to be in general unsolvable exactly [HI [14] . We also note that none of the known shape- invariant 
superpotentials listed, e.g., in refs. OUT], leads to exactly solvable partner Hamiltonians that can be experimentally 
implemented for molecules in nonresonant fields. For reference purposes, we provide a summary of the relationships 
among exactly solvable, shape invariant, and Infeld-Hull factorizable potentials in Appendix |A"| 

V. SUPERSYMMETRY AND SHAPE INVARIANCE OF THE COMBINED-FIELDS HAMILTONIAN 

IN THE FIELD-FREE AND STRONG-FIELD LIMITS 

In the field-free limit, w, Aw — > 0, the 3D Schrodinger equation reduces to one for a rigid rotor: 

JV(0) = Ei/>(6), (37) 
whose solutions are spherical harmonics, tp(9) = Yj m (9, 0) pertaining to eigenenergies Ej = J(J + 1). 
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FIG. 2: Eigenstates of a molecule in combined fields for different values of /3, with exactly solvable cases shown in blue. The 
equidistant energy levels in the strong-field limit, /3 — > oo, are labeled by the 3D librator quantum number n and are shown 
schematically. See text. 



(191, (20): 



We obtain the expressions for the superpotential and partner potentials of a free rotor by setting /? = in eqs. ( 18 1, 

W{9) = - (m + - J cot (9, (38) 



sin 2 (9 



m(m + 1) — - 



1 

4' 



and 



Vi m) (9) 



(m + l) 2 - i 



sin 



m(m + 1) 



1 

4' 



(39) 



(40) 



This is an obvious result as for J — m, the spectra of two rigid rotors with quantum numbers m and m + 1 coincide 
except for the ground state [12l [20] . 

Potentials (39) and (40) are seen to be shape-invariant, with a = m, ax = p(clq) = a + 1, and g(a n ) = a n (a n + 1), 



and the superpotentia 
eq. pi: 



of eq. (38) satisfies eq. (36). The eigenenergies of the Hamiltonian i/_ are obtained from 



E~ =n(n + 2m+l), (41) 

which for n = J — m coincides with the rigid rotor spectrum, shifted by —m(m + 1). The energy spectrum and 
wavefunctions of the supersymmetric partner potential V+ can be obtained via eqs. (23) and (24). 

In the strong-field limit, w, Au; — » oo, the angular motion of the molecules is confined near 9 = and potential ^ 
can be expanded in powers of 9. On retaining terms up to second order, the potential can be reduced to that of a 3D 



angular oscillator (librator). The corresponding Schrodinger equation ^ takes the form: 

f(9) = Ef(0), 



d 2 k 2 9 2 £{l + 1) 

de 2 + ~T~ + e 2 



where the eigenenergies are given by 



If 3 
E n = E + uj + Aw + - = k I 2n + £ + - 



n = 0,1,2, 



(42) 



(43) 



and k 2 = 2{uj + 2Aw), £ = (±m - 1/2). 

The analytic eigenfunctions of the angular harmonic oscillator can be expressed in terms of Laguerre polynomials [5] , 



f n (0) = N0 l+1 L^ +1 ^(x)exp(-x/2), x = \k9 2 



with 



N 



KX £/2+3/4 



(f) 



2n\ 



T{n + £ + 3/2) 



1/2 



(44) 



(45) 



and they vanish at x — > and decay exponentially at x — > oo 

e cons 

W{9) 



The corresponding superpotential can be constructed from the ground state wavefunction, eq. ( 44 ) , with n = 

f (6) k9 £+1 



MO) 2 

which leads to the supersymmetric partner potentials 

k 2 6 2 £{£+!) 



V-{6) 



n [£ + 



(46) 



(47) 



V+(0) = 



(£+!)(£ + 2) 

e 2 



The potential of eq. (47) coincides with the potential of eq. (42 ) shifted by the ground state energy, E = —n (£ + 3/2). 
Potentials (471 and (|48 ) are shape-invariant, with the parameters q,q = £, a% = p(clq) = clq + 1, and g(a n ) = 2na n 



(48) 



and superpotential (|46[) satisfies eq. ( 36 ) . The eigenenergies of the Hamiltonian H- can be obtained from eq. ( 34 1 : 



E n = 2nn, 



(49) 



which coincides with the eigenenergies given by eq. (43) shifted by — k£ — 3/2. The eigenenergies and eigenfunctions 
of the superpartner potentials also fulfill eqs. ([23]) ( 25 1 and represent a complete analytic solution to the problem. 

Figure [2] shows the energy levels of a molecule in combined fields for different values of the field-strength parameter 
j3. In the weak-field limit, /? —¥ 0, the energy levels approach those of a free-rotor, which is solvable exactly for all the 
eigenstates. For nonzero but weak fields, /3 = 1, the levels become split in J + 1 components due to the Stark effect. 
In this case, the SUSY partner Hamiltonians are not shape-invariant, and the problem is analytically solvable only for 
the "stretched states," corresponding to J = m. With increasing interaction strength, the energies of the stretched 
states come closer to one another and, in the strong-field limit, j3 — > 00, coalesce into the ground state level of the 
3D harmonic librator. In the strong-field limit, the supersymmetric problem becomes shape-invariant again, and is 
exactly solvable for all eigenstates in closed form; the equidistant levels are infinitely degenerate and separated by an 
energy difference of (2w + AAuj) 1 ^ 2 . 



VI. APPLICATIONS 

The analytic wavefunctions of molecules in combined fields make it possible to obtain their properties analytically 
as well. The space fixed dipole moment, fiz, is given by the orientation cosine, (cos 9) = (ip(6)\ cos 8\i/j(8)), and for 
the exact wavefunction of eq. (31 ) can be evaluated in closed form, 



Hz/n = (cos 6) = - 

J- 



m+1 



(2/3) 



h3/2 



(50) 
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where I n (z) is a modified Bessel function of the first kind [T3]. Fig.[3]shows spaced fixed dipole moments corresponding 
to the states |J = m,m;/3) for several values of m as a function of the j3 parameter. The value of [iz rapidly 
increases with (3. For instance, for m = 0, it rises from only 0.54/^ at f3 = 1 to 0.83/z at (3 = 3. In the case of the 
much studied 40 K 87 Rb molecule, which possesses a dipole moment fi = 0.589 Debye and a polarizability anisotropy 
Aa = 54.21 A 3 [HJ relatively weak fields of e = 38 kV/cm and / = 1.75 ■ 10 9 W/cm 2 (corresponding to ft = 5) 
give rise to a strongly oriented ground state with [iz = 0.9/i. This value of (cos 8) corresponds to the molecular axis 
confined to librate within ±26° about the common direction of the fields. 

The alignment cosine, (cos 2 8) = (if>(0)\ cos 2 8\ip(8)), characterizes the molecular alignment along the Z axis. For 
the wavefunction of eq. (31 1, the alignment cosine takes the analytic form, 



(cos 2 8) = ^Q^i O m + 7 A P 2 ) + 0*1 (i m + 5/2; (3 2 ) 



2 Q F 1 (;m + 3/2;l3 2 ) 



(51) 



with Qpi(;a;z) — oFi(;a:z)/T(a) a regularized confluent hypergeometric function [T3]. Fig. [3] shows (cos 2 8) of the 
states I J = m, m; (3) for several values of m as a function of the j3 parameter. The |0, 0; (3) state exhibits quite a strong 
alignment with the alignment cosine rapidly approaching with increasing f3 the value of 0.8, which corresponds to a 
libration of the molecular axis about the polarization vector of the radiative field with an angular amplitude of 27°. 
The expectation value of the angular momentum is related to the orientation cosine, eq. (50 1, via 



hn+l 



(2/3) 
/2 (2/3) 



/3 (cos 8), 



(52) 



We note that the dependence of (J 2 ), shown in Fig. [5J becomes asymptotically linear in j3 for all the values \m\, cf. 
eq. ((52J. 

By making use of eqs. ^ and ((FT}, one can show that the following condition for the expectation values is satisfied: 



sin 



2(3 (to + l)(cos0) - (3 2 (cos 2 8) = E , 



(53) 



with the energy Eq given by eq. (16 1. 

Table [I] lists, for reference purposes, the above analytic expressions for the space- fixed dipole moment, alignment 
cosine, and expectation value of the angular momentum for different values of \m\. 



TABLE I: Analytic expressions for the space-fixed dipole moment fiz/n = (cos 6), alignment cosine (cos 2 ( 
value of the angular momentum (J 2 ) obtained from eqs. (501, (511, and (52 1 for different values of \m\. 



and expectation 



■0) 



I coth(20) 
2/3 2 
3 8/3 
& 1 tanh(2/3) - 2y3 

15 8/3 2 



2/3 3 
2/3 coth(2/3) - 1 20 



3 (4/3 2 - 6/3 coth(2/3) + 3) 3 



3 + 4£ 2 - 6/3 coth(2/3) 



/3coth(2,3) - 



2/3 2 
2/3 coth(2/3) ■ 



3(3 + 4/3 2 - 6/3 coth(20)) 



(16/3 4 + 180/3^ 4- 105) - 10/3(8/3^ + 21) coth(20) 2/3(4/3 4 4- 95/3^ + 210) coth(20) - (56/3 11 4- 3750^ 4- 210) -4/3(7/3 z 4- 15) coth(20) 4- 2(4/3 4 + lift* 4- 15) 



4/3 2 (4/3 2 4- 15) coth(2/3) - 6/3(8/3 2 + 5) 



2/3 J (4fl2 + 15 ) coth (2/3) - 30 2 (80 2 4- 5) 



23(4/3 2 + 15) coth(2^) - 3(80 2 4- 5) 



VII. CONCLUSIONS 



By invoking supersymmetry, we found a condit ion under which the molecular Stark effect problem becomes exactly 
solvable. The condition, Aw = 4 („" +1 )2 , cf. eq. (17 1, connects values of the parameters u) and Aw that characterize 
the interaction strengths of a polar and polarizable molecule with collinear electrostatic and nonresonant radiative 
fields. The exact solutions are obtained for the |J = to,to;cj, Aui) family of "stretched" states. 

We also considered the field-free and strong-field limits of the combined-fields problem and found that both exhibit 
supersymmetry and shape-invariance, which is indeed the reason why they are analytically solvable. 
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FIG. 3: Ground-state energies Eo (in units of the rotational constant B), space fixed dipole moments \iz / M = {cos 6), alignment 
cosines {cos 2 8), and expectation values of the angular momentum (J 2 ) for different |J = m,m;/3) states as a function of the 
interaction parameter f3. 



By making use of the analytic form of the | J = m, m;cj, Auj) wavef unctions, we derived simple analytic formulae for 
the expectation values of the space-fixed electric dipole moment, the alignment cosine, and the angular momentum. 
These key characteristics of molecules in fields are summarized in Table [I] We also derived a "sum rule" which yields a 
formula for the corresponding eigenenergy, eq. (53 1, in terms of the above expectation values. The analytic expressions 



obtained open a direct route to engineering molecular states with preordained characteristics. 

Interestingly, it is possible to glean the reason as to why the exact solution of eq. (pi is obtained for only one 
particular relation between the field strength parameters, eq. (17), from the semiclassical (WKB) approximation |23| . 



The eigenfunction of a ID Schrodinger equation assumes the WKB form f(9) oc exp[iS(9)], with S(9) the action of the 
underlying classical system. For the ground state of the potential ( fl9| , the action satisfies the differential equation, 
S'{9) 2 — iS"(9) — Fi 1D , whose solutions are obtained by expanding S(9) in powers of K. It turns out that in the case 
of the combined field strengths connected via eq. (17), the series converges to the following exact expression, 



S(9) = — [2/3 cos 6 + (2m + 1) ln(sin t 
2i 



(54) 



which, when substituted into /(#), yields the exact, closed form wavefunction, eq. (27). 

We note that the exact ground-state wavefunction |J = 0,m = 0;uj — 0, Aw) can be also obtained as a "curious 
eigenproperty" by the method outlined by von Neumann and Wigner in 1929. They showed that by imposing the 
integrability condition on the sought wavefunction, a class of potentials could be derived that support a localized 
bound state embedded in the continuum P4lt26j . 

An extension of the method described in Section |III| might furnish other types of exact wavefunctions for the case 
of more than two combined fields and for non-collinear field geometries. 
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Appendix A: Exactly solvable, shape invariant, and Infeld-Hull factorizable potentials 

Supersymmetry, shape invariance, exact solvability, and the factorization method are often studied concurrently in 
the literature. In this appendix we briefly survey the literature studies on the relationships among them. Figure [4] 
summarizes the results. 

Every one- dimensional potential with a ground state in close form admits SUSY. Given a one dimensional 
quantum system with at least one bound state, one can find a partner Hamiltonian which has exactly the same 
discrete spectrum except for the ground state energy of if _ [37J 12E] • 

Every Infeld and Hull factorizable potential is shape invariant but the converse is not true. Shape invariance 
offers, in general, more than the factorization method, since the factorization method treats only the translational 
shape invariance j28, '2S]. The table prepared by Infeld and Hull 12 , although complete for most purposes, is not 
the most general table possible, as it does not include the most general solution of the Ricatti equations, worked out 
in Ref. [29]. 

Shape invariance and normalizability are sufficient but not necessary conditions for exact solvability. Gendcn- 
shtcin suggested that all exactly solvable potentials must be shape invariant [30] . but many counter examples to 
this conjecture were later constructed, for instance the Natanzon class of potentials which are, in general, not shape 
invariant [TH1 133 122] • We also note that in order to be exactly solvable, the shape-invariant potentials should admit 
normalizable solutions. 




Sup eisymme trie 
Exactly Solvable 
Shape Invariant 
Infeld and Hull Factorizable 

FIG. 4: A Venn diagram showing the relationhips among supersymmetric, exactly solvable, shape invariant, and Infeld & Hull 
factorizable potentials. 
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